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Abstract 

We investigate violations of the fluctuation-dissipation theorem in two classes of trap 
models by studying the influence of the perturbing field on the transition rates. We show 
that for perturbed rates depending upon the value of the observable at the arrival trap, 
a limiting value of the fluctuation-dissipation ratio does exist. However, the mechanism 
behind the emergence of this value is different in both classes of models. In particular, for 
an entropically governed dynamics (where the perturbing field shifts the relative population 
of traps according to the value of the observable) perturbed rates are argued to take a 
form that guarantees the existence of a limiting value for the effective temperature, utterly 
related to the exponential character of the distribution of trap energies. Fluctuation- 
dissipation (FD) plots reproduce some of the patterns found in a broad class of glassy 
systems, reinforcing the idea that structural glasses self-generate a dynamical measure 
that is captured by phenomenological trap models. 

1 Introduction 

Glasses are widespread materials that yet resist a basic comprehension from many points of 
view. Although several aspects concerning their macroscopic behavior have been modeled, 
showing qualitative agreement with experimental results, still our understanding of the 
relevant physics at the mesoscopic scale remains largely speculative. Long timescales (i.e. 
several orders of magnitude larger than the microscopic timescales) associated to physical 
processes occurring at very small length scales (of the order of several nanometers) are 
thought to be common. The main difficulty behind a theoretical description of glassy 
dynamics is caused by the great disparity between timescales and length scales, where 
neither hydrodynamic approaches or standard nuclcation theories can be applied. It is 
for this reason that trying to understand the dynamical behavior of glasses remains a 
theoretical challenge. Its comprehension may open new routes to the understanding of 
a large class of new physical phenomena, where such a strong disparity of temporal and 
spatial scales is relevant. 

The glass transition problem has been studied from many perspectives. In the ab- 
sence of a well established basic microscopic theory modeling offers many possibilities. A 
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successful category of models for the glass transition are trap models They have 

the advantage of capturing some of the relevant features of glassy relaxation within a 
tractable formalism. Trap models have been widely investigated (for a review, see 0) 
and found to reproduce pretty well several experimental features such as the decay of the 
thermoremanent magnetization in spin glasses IKS| or the super-activated behavior of 
the viscosity in glasses jS] . Glassy systems are characterized by aging and violations of the 
fluctuation-dissipation theorem (FDT) (for recent reviews see [E1EJ). FDT violations can 
be quantified by defining the fluctuation-dissipation ratio (FDR), 

ds 

where T is the temperature of the bath (the Boltzmann constant ks is set to 1) and 
C(t, s), R(t,s) are the two-times correlation and response functions. The FDR allows to 
define also the effective temperature, 

T eS (t,a)=T/X(t,a) . (2) 

In equilibrium FDT is satisfied so X = 1 and T e g(t, s) = T. These definitions would be 
devoid of any interest if it were not for the fact that a broad category of glassy systems 
(those exhibiting a one-step replica-symmetry breaking transition in their corresponding 
mean-field theory) show, in their long-time regime, that X(t,s) — > X(s), i.e. the FDR 
converges to a non-trivial function of the lowest time s. More precisely, these systems 
show the existence of two-timescales depending upon the value of the ratio (t — s) / s in 
the asymptotic regime t,s >> 1. If (t — s)/s << 1 then X(t,s) — 1 and the systems is 
locally equilibrated (short-time sector). However, if (t — s)/s ~ 0(1) then X(t, s) —> X(s) 
(long-time sector) being independent on the precise value of the ratio t/s (provided it 
is finite). Hereafter, we will refer to this category of models as one-step (IS) models as 
they show a single jump of the value of X(t,s) from 1 to a different value (generally 
smaller than 1). The existence of X(s) shows that the aging state at time s can be 
characterized by a non-equilibrium parameter that could be identified as the fingerprint of 
a non-equilibrium temperature, T e fj(s) = T/X(s). From a theoretical point of view, the 
conditions necessary for the existence of X(s) in the long-time sector have not been found 
in full generality. These conditions could be quite strong, though, and the possibility has 
been raised that only the limiting value X oc (s) = lim£_>oo X(t, s) (rather than X(s)) is 
physically meaningful . Most of the numerical studies done on structural glass models 
seem to be in qualitative agreement with the IS scenario pj. 

Usually, FD violations are schematically depicted in the form of FD plots where the 
integrated-response function x(^ s ) — / R(t,t')dt' is plotted as a function of C(t,s). For 
IS models FD plots are made out of two linear segments, one of slope —1/T covering the 
short-time sector, the other with slope —X(s)/T in the long-time sector. 

Recently, a debate has emerged whether trap models can reproduce the pattern char- 
acteristic of FD violations in IS systems. In particular, because dynamics within traps is 
usually not incorporated into the models, they should only reproduce the long-time sector 
where X(s) describes relaxational jumps among traps. Yet, the most recent results on this 
issue ^T] seem to disprove the existence of such IS features. 

The goal of this paper is to clarify this issue by showing that FD plots in trap models 
strongly depend upon the functional dependence of the transition rates on the field. Actu- 
ally, perfectly straight IS FD plots can be obtained provided the escape rate from a trap 

The statement has been done in the context of ferromagnetic models at criticality where the universal 
FDR is Xoo = lim MOO Jfoo(s). Several analytical works support the conjecture |§1 110|. The conjecture can 
be naturally extended beyond criticality to generic glassy systems by proposing that X QK1 (s) is waiting-time 
dependent. 
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is unchanged (within the linear-response regime) by the action of a external perturbation. 
The study of two different classes of trap models shows the existence of a value (s) for 
the FDR (JJJ in the limit t — > oo, yet the mechanism behind the emergence of this limiting 
value is different. For a dynamics governed by entropic effects there exists an effective 
temperature T ao (s) that remains finite in the limit T — > and is uniquely determined by 
the width of the exponential distribution of trap energies. This shows that f S behavior 
with X oa (s) proportional to T at low T is the fingerprint of a slow dynamics governed by 
entropic rather than thermally activated effects. 

2 Trap models: some definitions and useful expres- 
sions 

Following pQ we consider N independent traps labeled by the letter a, each trap character- 
ized by an energy E a and a generic observable that we will denote by M a . We consider a 
Markov dynamics among traps. To set the formalism as clear as possible we will consider 
time as discrete and define a dynamical history by the temporal series {a(t)} where t is an 
integer (and consecutive time steps are separated by At = 1). The dynamical evolution 
is specified by the probability P t (a) that the system stays in trap a at time t, and by 
the conditional probability W(a'|o;) that the system goes from trap a to a' in a single 
time step. We assume conditional probabilities to be time- independent. Eventually we 
will be interested in the continuous-time limit At — ► where conditional probabilities for 
a ^ a' are vanishingly small W(a'|a) — > W(a'\a)At, W(a'\a) being the transition rates. 
In what follows we will coin the term rate for W although properly speaking this term 
should be only adopted for its continuous-time analogue W. A master equation governs 
the dynamical evolution of the system, 

P t (a)=J2P t ^(a')W(a\a') . (3) 

a' 

The W will consist of two terms, the term W(a'|a) = A(a'\a) for a ^ a 1 describing the 
jumping rate from a to a' and the term W(a|a) = B(a) describing the rate of staying in 
trap a. Therefore we can write, 

W(a'\a) = A(a'\a)(l - 6 a ,a>) + B{a)6 a , a , (4) 

Transition rates are normalized J2 a W(a\a') = 1 and satisfy, 

A(a'\a) + B(a) = 1, Va . (5) 

a' 

The transition rates allow to define the escape rate W esc (a) from a given trap a, 

W csc (a) = J2 Ma'\a) = I - B{a) (6) 

Dynamics in trap models is specified once the rates W are given. For the time being we 
concentrate only in the zero-field case where dynamics is only governed by the energies 
of the traps. We consider two representative cases: the Bouchaud model (BM) 1 where 
the rate A(a'|a) only depends on the energy E a of the departure trap a, and the Barrat- 
Mezard model (BMM) 2 where the rate A(a'\a) depends on both energies of the departure 
and arrival traps, E a and E a > . Both models consider an spectrum of traps with energies 
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E a < and arc defined by their respective rates, 



i BM (a»^exp(R) (7) 

(BMM/„/u ^0 1 

A {a |a) = Nl + cMP(E a ,~E a) ) (8) 

where (3=1 fT and loq is a microscopic frequency that we will take inversely proportional 
to the discretization time At. The choices (|7I8|) are justified by the requirement that 
the detailed balance property is satisfied. The corresponding expressions for B(a) are 
determined by the normalization condition J^J ■ The energy spectrum of the ensemble of 
traps is defined by, 

N 

X 

a=i 

Both models consider an spectrum of states with exponentially distributed energies, 



i ,y 

p(E)=kY, S ( E - E °) ( 9 ) 



p(£) = i_exp(^) E<0 . (10) 
1 g yl g / 

There is an important difference between the dynamics in the BM and the BMM in the low 
T limit. While dynamics is completely arrested in the BM (the escape rate (JHJ) vanishes 
exponentially fast with 1/T) it does not in the BMM where the escape rate is finite even 
at T = 2 . This makes the BM model more suitable to describe glassy relaxation in a 
landscape of traps separated by energy barriers while the BMM is better suited to describe 
relaxation over entropy barriers. 

In the presence of an external field h coupled to the observable M a the rates (|7I8|I get 
modified. The simplest possibility is to assume that the functional form of the rates (|TIHf> 
remains unchanged while replacing 

E a ~^E a - hM a . (11) 

We will refer to this choice as the purely activated rates. A priori, however, other choices 
maybe equally valid. For instance, let us consider the following rates in the presence of a 
field 3 , 

A h (a'\a) = A(a'|a)ex P (-/3/iO«M Q - jM a ,)j (12) 

where [i and 7 are arbitrary parameters. In this form these rates generally violate detailed 
balance. Although this choice may seem unjustified, violation of detailed balance may 
be a necessary ingredient for trap models to reproduce non-equilibrium aspects of glassy 
systems 4 . 

The rates (JT2J were considered in 0] in the case where 7 = 1 — fx and detailed balance 
is preserved. However, the implications and physical meaning of the term 7-M Q ' in (|12JI 



2 This difference between both models appears explicitly in the time decay of the energy in the glassy phase 
T < T g where the Gibbs energy distribution ceases to be normalizable. While in the BM the energy decays 
as E(t) ~ — Tlog(t) (therefore relaxation arrests at T = 0), in the BBM it decays like E(t) ~ — T s log(t) for 
T — > being T independent in that limit. 

3 From now we will use the subscript h to refer to field-dependent rates. For zero-field rates we will drop off 
the subscript h. 

4 Actually, trap models can be justified only as a coarse-grained description where the original phase space 
is partitioned into components or regions, each of these regions corresponding to a trap (for a discussion 
of this aspect see Sees. 3.3 and 4.2 in 6 ). However, at difference with the original microscopic dynamics, 
the projected dynamics among traps can violate detailed balance (a property that must be satisfied by the 
microscopic dynamics for the system to equilibrate). 
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went unnoticed. Indeed, the choice of purely activated rates (p, = 1,7 = 0) has occupied 
the majority of studies of FDT violations in trap models @J^2E|- Here we go beyond 
and show that the existence of the term ^M a i is crucial for the emergence of IS behavior in 
FD plots. We will give also a physical interpretation of this term in a relaxational scenario 
governed by activated and entropic barriers. Our study will consider both the BM and 
the BMM. 

Before delving into each of the two class of models we present some useful formulae 
needed to derive the FDR. Some of these computations have been addressed in previous 
works 01 1111 118) , here we reproduce them in more detail. We are interested in comput- 
ing the response function R(t, s) and express it in terms of derivatives of the correlation 
function C(t,s). They are defined by, 

C(t,s)=<M(t)M(s)> ■ R(t,s) = ± S< s A h I ^ > 6(t-s) , (13) 

i.e. the response function is the change of the expectation value < M{t) > when a impulse 
field coupled to M is applied at a previous time s. For simplicity we will assume that 
C(t,t) is normalized to unity. If Wh{a'\a) stands for the rate in the presence of a field 
then we have, 



R(t,s) = — Ps{a s )M at J] W{a r+1 \a r )A h W(a s+1 \a s ) (14) 

where 



At 

a>t,...a. s r— s+1 



A h W(a'\a) = lun W ^- W ^ (15) 

h^O h 

Using (1 1I5[) we can rewrite (|T3|l as follows, 

A h W(a'\a) — a(a'\a) - j^y^ a{a" \a) 5 

a', a. 

(16) 

a" 

with the definition, 

( dA h {a'\a) \ 
a(a|a) = l dh A=o ' 

The next sections describe how to compute a(a'\a) in both the BM and BMM and how 
an expression for R(t, s) in (|14|l can be obtained. 

Another important concept related to the measured observables M a is neutrality |lll!5|. 
The probability density of being at time f in a trap with a value M of the observable, 
conditioned by its energy to be equal to E, is given by, 

2^ a Pt(a)d{E - E a ) 

M is said to be neutral if, assuming that P t= o(M\E) — P t= o(—M\E), this property stays 
valid at all times Pt{M\E) = P t (—M\E). As stated, this condition of neutrality is quite 
stringent. It can be relaxed provided some conditions are met. For sake of simplicity and 
to successfully convey the main message of this paper, we will assume that the observable 
density g{M) is symmetric 5 , 

1 - 

9i M ) = Jj E 5 ( M - M «) = 9(~M) (19) 



3 A common choice is a Gaussian distribution of zero mean and unit variance such that C(t,t) = 1 holds. 
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Eq. (|T§|> together with the initial condition P t =o(M\E) = P t=Q (— M\E) is enough to 
guarantee neutrality. Neutrality allows to derive a simple relation for the lowest-time 
derivative of the correlation function C(t, s), 



dC(t,s) 



ds At 

a t ,..,a 



1 t_1 

- Ps(a s )M at l[W(a r+ i\a r )(M as+1 - M a J (20) 



If we replace (@J for ^(c^+i la^) in l|20|l then neutrality implies, 

Y / P s (as)A(a s )M as =0 (21) 

because, for a fixed value of the energy, an identical number of traps contribute to (|2 1 1) 
with positive and negative values of M . Therefore, 



dC(t,s) 



1 

,-> s . 'At ^ P s(u s )M at J] W(a r+1 \a r )A(a s )M as+1 . (22) 

Qt , .. ,ot s r— s-\-l 

3 The Bouchaud model (BM) 

In the BM the unperturbed rates A BM (a'\a) are given by and depend only on the 
energy of the departure trap so we will adopt the notation A BM (a'|a) = ^4 BM (a). Let us 
consider the rates lfl2"|l . 

Af M (a'\a) = A BM (a)exp(-/3/i(/iAf Q - -yM a ,j) (23) 

It is easy to verify that in this case a{a'\a) is given by, 

a BM {a'\a) = —[3A BM (a) ([J.M a - 1 M a ,) (24) 

where [i and 7 are two constants which we will assume can take any value (positive or 
negative). The parameter 7 appearing in the rates Q24J) only depends upon the value of 
the observable M at the arrival trap. At linear order in h, the change of the escape rate 
W csc (a) © is given by 

W^ sc (a)-T^ csc (a) ~ -vPhNA BM {a)M a + 0(h 2 ) = -uj nf3hM a exp{pE a ) + O(h 2 ) (25) 

and does not depend on 7. Therefore, while the parameter /i directly modifies the lifetime 
of a trap, the parameter 7 does not modify it (up to linear order in the perturbation). 
Simply put, while the parameter fi directly affects the height of the energetic barrier to be 
surmounted (through thermal activation) to escape from a trap, the parameter 7 selectively 
enhances the rate toward the destination trap through a different activated mechanism. 
In Sec. [3 we present a physically motivated explanation for this dependence. 

Note that only if 7 = 1 — /i detailed balance is preserved in l|24() . Inserting (|24|l in 
fTCim we get, 

A h W(a'\a) = -/3A BM (a) [( P M a - 1 M a ,) - N(fiM a - 7 M)<5 Q ^] (26) 

where M — (1/N) J2 a M a is the average magnetization of traps. The following mathemat- 
ical steps are followed. We insert (|26|l in (|14|) and use the identity (@J . If the observable 
M a is neutral then from H19(l we have M — 0. Using l|22(l we get, 

r dC(t,s) dC{t,s) ^ 
fl(t,.) = -0[ M _ gr -- 7 _ g _| (27) 
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where the continuous limit At — ► has been taken 6 . Replacing in (JTJ we get for the 
FDR, 

/ 9C(t,s) \ 

Z(t, S )= 7 - / x(^ T j (28) 

Several comments are in order. First, the present analysis covers regimes where detailed 
balance is violated. We already justified this choice (see footnote 0J. In fact, in a time- 
translational invariant regime, X = 7 + /1 which is generally different from 1. Only if 
detailed balance holds (7 = 1 — fi) then X = 1. Second, X is finite in the T — > limit. In 
other words, if a physically meaningful effective temperature exists, then it has to be 
proportional to T and vanish at T — 0. Let us consider the following aging form of C(t, s), 

poo .j. roc 

C(t,s)=j d\w(X)C\ [—^-j J d\w(X) = l (29) 

where C\(0) = 1, C\(oo) — 0. Provided w(X — > 0) — > then it can be easily proved that 
for t ^> s, 

X(M)^7-e>(|) (30) 

and therefore 

Xoo = lim X(t,s) =7 (31) 

i — >oo 

i.e. the FDR asymptotically converges to a non-trivial quantity (7). The main conclusion 
of this simple calculation is that for rates given by Ij24(l there is a universal dependence 
of the FDR (in its long-time limit) on the parameter 7 describing how transition rates are 
perturbed at the linear response level (i.e up to linear order in h). Note that for \x = 0, 
X(t, s) = 7 (|28p. the FD plots are full straight lines as in the IS scenario and the escape 
rates are unmodified at the linear order in the field (|25() . 

The validity of (|30(1 relies on the fact that the system is aging. Were the system in a 
stationary state then trivially X(t,s) = 7 + /1 for all times t, s. Therefore, for stationary 
non-equilibrium systems X x = 7 + fi and an additional dependence of the limiting value 
Xqo on /i is expected. 

In general, for models with full aging C(t,s) = C(t/s), X(t,s) in lt2"5|) is a single 
function of C yielding limiting FD plots For instance, the BM is characterized by 

a non-normalizable distribution for T < T g where the correlation function C(t, s) shows 
full aging. In Fig. ^ we show FD plots for the BM for a given value of 7 and different 
values of \i. Note that for ijl = a IS FD plot is obtained characteristic of systems with a 
single timescale. The purely activated case \x — 1,7 = gives X^ = as has been found 
inEl. 

4 The Barrat-Mezard model (BMM) 

As we already said in Sec. El after (|10|l the main difference between the BM and the 
BMM is that the escape rates in the later are temperature independent in the low T limit, 
therefore slow dynamics is governed by entropy barriers. Therefore, the BMM appears as 
the phcnomenological description of the entropic relaxation of several microscopic models 
such as the Backgammon model [T3] or the oscillator model |15| . 

In this section we want to show that similar conclusions (|25I31|I can also be inferred 
from the BMM JHJ in the low-temperature limit T — > where a purely entropic mechanism 

6 Equation 127H was first derived in |1] but for the case where detailed balance holds. Note also that the 
symmetry condition there expressed (our i!9l ) is not enough for the validity of 1271 as full neutrality as 
expressed by the supplemental condition P t= o(M\E) — P t=0 (— M\E) must hold. 
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Figure 1: FD plots for the BM at T = 0.5 (T g = 1 with g{M) a Gaussian distribution of unit 
variance) at 7 = 1 and \x = 0, 0.5, 1, 2, 3 (from bottom to top). Note that only for /i = there 
is a full straight line, yet in all cases the limiting slope when C — > is the same, = 7=1. 
It can be proved that the value of x(C) i n the limit C — > is given by ^7 + (i,/T g . 

of relaxation is at work (in contrast to the BM case where escape rates are thermally 
activated). The analysis of the FDR in the BMM goes along the same lines as in the 
preceding section. We start from the unperturbed T = rates, 

A BMM (a'\a) = ^6(E a -E a ,) . (32) 

We remark two aspects of these rates: 1) As these are zero-temperature rates they do not 
satisfy detailed balance and 2) the rates are independent of the energy of the departure 
and arrival trap provided that the energy decreases. Nevertheless the escape rate © from 
a given trap will depend on its energy through the relation, 

W csc (a) = ^J2 9 ( E " - E <*')= ^ / " P{E')dE' (33) 

J-00 

and will decrease as the energy of the trap E a decreases. Now we need to specify how the 
rates W change in the presence of a field. Inspired by the corresponding expression for 
the BM H23|l we will assume the following factorized form for the perturbed rates in the 
BMM 

Af MM (a'\a) = A BMM (a'\a) exp(-/i(/zM a - 7^)) . (34) 

where A BMM (a'\a) is given by (J22J. This functional dependence assumes that transitions 
to new traps are enhanced towards larger values of the observable M a i , however arrival 
traps must have energies below that of the departure trap E a i < E a . This form of the 
perturbed rates does not correspond to the case of purely activated rates where l|3"2"|) are 
modified according to the shift In fact, were this the case then there would be 

correlations between energies and observable values, and transitions with E a > > E a would 
occur depending on the difference M a i — M a . The factorized form of the rates l|34|l as 
well as the enhancement of rates towards traps with larger values of M a are not easy to 
justify a priori. They have to be seen as an approximate solution to the real perturbed 
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rates that takes into account the dependence of the rates upon the observable value taken 
at the arrival trap 7 . 

The enhancement of rates towards traps with larger M a i can be explained as follows. 
Let us assume that p(E) is a monotonically increasing function of E (e.g. as in I|1U|)1. The 
main effect of the field on a shell of traps that have a given energy is that the population of 
traps with M > tends to increase relative to those with M < 0. The shift of populations 
can be computed as follows. Given the unperturbed density p^ 2 \E,M) = p(E)g(M) we 
can write, 

p { *\E,M)= dE' p(E')g{M)5{E - E' - Mh) = p{E + Mh)g(M) (35) 



= p(E)g(M)(l+ dl0e ^ E)) Mh + 0(h 2 )) =pW(E,M)exp( 7 (E)Mh) (36) 

the final expression being valid up to linear order in h and where we have introduced the 
function j(E), 

7(£) = W>>=^ , ( 37> 

where S(E) = log(p(E)) is an entropy an 7 has the dimensions of the inverse of a temper- 
ature. Analogously, if we define an entropy S {2) {E,M) = \og(pW(E,M)) then we have, 
up to linear order in h, S^ } (E, M) = S (2) {E, M) + j(E)Mh. The schematic behavior of 
{E, M) is depicted in Fig. For an exponential distribution of states (fTOf) we have, 

7(E) ^ 7 = ^ (38) 

1 9 

Although the density Ph(E, M) changes at linear order in h, the symmetry property (|19J) 
ensures that the total density Ph(E, M)dM only changes at the next order h? . We 
stress that the parameter 7 (expressing the dependence of the perturbed rates upon the 
value of the observable evaluated at the arrival trap), rather than p, is the only term 
describing the entropic unbalance of traps. From Ij6l34|l we obtain for the escape rates, 
at linear order in the field, 

WT(«)^«P(^) (39) 

where we used the result M = and l|33l) . Calculations identical to those shown for the 
BM in Sec. El lead to, 

dC{t,s) dC(t,s)- 



R(t,s) 



V — HI 7- 



(40) 



dt ' ds 

Note that, contrarily to l|27|l . the response is now finite in the T = limit, therefore 
the FDR vanishes linearly with T as T — > 0. Rather than X it is useful to consider the 
effective temperature P|l. 



T e s(t, s) 



7 - M 



dC(t,s) 

dt 
dC(t,s) 

ds 



(41) 



7 A contribution to the perturbed rates of the form 134H seems to us essential for the emergence of non-trivial 
effective temperatures. After completion of this work, calculations by Sollich (see note added at the end of the 
paper) have shown that purely activated rates are equivalent to l|34|l with p = 2/T 9 ,7 = where there is no 
dependence of the modified rates upon the observable value at the arrival trap. 
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Figure 2: Illustrative picture for (E,M) (Gaussian distribution g(M) of unit variance and 
T g = 1) as a function of M for a fixed value of the energy, h = (continuous line) and h = 0.1 
(dashed line). The difference between the two curves is linear with h and the slope is equal 
to h/Tg = 0.1. 



Again, under the same assumptions (|29|l and using l|38|l we get in the limit t — > oo, 

Too (a) = lim Teff(t, s) - - = T g (42) 

which is independent of s and Therefore, in the limit t — > oo the effective temperature 
coincides with T s which is the width of the exponential distribution of states through 
(j37l38|l . The change of the escape rate (© induced by the field is given by 

W c h sc (a) - W csc (a) = -cu of ihM a exp(^) + 0(h 2 ) (43) 

and the escape rate, as for the BM, only depends on the value of /i. Again, full straight 
FD plots are obtained for /i = when the escape rates (or trapping times) are unchanged 
by the field. For the BMM, computations of I|41I42|) can be made explicit due to the 
simplicity of the expression for the correlation function, C(t, s) — s/t for (,s>l 2 j. For 
the integrated response function we get, 

X (t, s) ee x(C) = f dt'R(t, = 7 (1 - C) + |(1 - C 2 ) (44) 



and the effective temperature, 

1 = d X {C) 

Teff(C) dC 



= 7 + nC (45) 



giving for ji = straight FD plots characteristic of IS behavior, x(C) = 7(1 — C). Note 
that, in the \x = case, x(C = 0) = 7 = 1/T g which is the value of the susceptibility at 
the transition temperature. Fig. shows some typical FD plots. 



4.1 Numerical test of the BMM solution with rates (|34j) . 

To verify the correctness of the previous solution we have done a numerical simulation of 
the BMM with rates (|34J) . We have considered a Gaussian distribution of observables 
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Figure 3: FD plots for the BBM at T = (T g = 1 with g{M) a Gaussian distribution of 
unit variance) at 7 = 1 and fi = 0,0.5,1,2,3 (from bottom to top). Note that only for 
(i = there is a full straight line, yet in all cases the limiting slope when C — > is the same, 
Too = T g = 7 = 1. 

g(M) = (27r) -1 / 2 exp(-M 2 /2) and T g = 1. Let E, M be the energy and trap at the 
departure state. The procedure is as follows: 1) a trap with observable M' is chosen 
with a probability proportional to exp(M'h/T g )g(M'), i.e. M' — r + h/T g where r is 
a Gaussian (zero mean and unit variance) random number; 2) the lifetime of the trap, 
according to (|39[) . is given by r = (1/wo) exp(— E + Mh); 3) The arrival energy is given 
by E' = E + log(r') where r' is a random number uniformly distributed between and 1. 
We start from the disordered initial state E = 0, M = r where r is a Gaussian random 
number as before. We then run two independent simulations by iterating many times the 
previous steps (one run is at h = for all times, the other run is at h = until time s and 
later a small field h is switched on) and generate correlations and susceptibilities (|44|) at 
different times t, s. For t, s > 1/u)q the asymptotic aging solution is obtained as well as the 
limiting FD plots Ij44(l . We verified the correctness of the solution C(t,s) at zero field. 
The numerical results are shown in Fig. 0] (right) together with the analytical prediction 
<|38|) with 7 = /j, = 1 (continuous line). 

5 Discussion 

The entropic mechanism described in the preceding section for the BMM provides an expla- 
nation on the emergence of an effective temperature given by the width of the exponential 
density of trap energies 8 . However, is it possible to extend that explanation to the BM 
where such entropic mechanism appears difficult to reconcile with the explicit activated 
behavior of the rates (|23|l ? In fact, there is an important difference between the BM and 
the BMM. While in the later the perturbed rates are not thermally activated (see (13411 ) 
leading to a value of Xoo proportional to T at low T, in the former the perturbed rates 
are thermally activated 12- il) leading to a value of that is temperature independent 
at low T. 

The following explanation, albeit speculative, might be useful. Consider a glass where 
relaxation is driven by the nucleation of cooperatively rearranging regions. The nucleation 

8 In the mean-fieid theory of spin giasses such origin of the effective temperature has been aheady re- 
marked 117) . Yet in these modeis activated processes are compieteiy suppressed. 
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Figure 4: Numerical test of (jMJ). FD pl° ts f° r the BBM with 7 = = T ff = 1, u = 1 and 
two waiting times s = 10 (circles) and s = 1000 (crosses) in the asymptotic regime s >> 1/uq. 
The continuous line is the analytical result (|44j) . 

of such a region requires the simultaneous activation of M independent processes, each 
one characterized by a thermally activated rate w 9 . In the presence of a perturbation of 
intensity h the rates get modified Wh ~ wexp(—/3h(ji), where w is the rate at zero field 
and fi is a microscopic parameter. The modified rate for the nucleation of the region is 
wff ~ w M exp(— (3h/J,M). If a, a' denote the configurations of the system before nucleating 
and after nucleating respectively, then M can be identified with the net change of the 
observable M a — M a i after nucleation takes place. The overall rate for that nucleating 
process is exp(—(3hfJ,(M a > — M a )) which coincides with (|23|1 with 7 = /i. 

This scenario could describe diffusive processes (such as low T coarsening) where partial 
equilibration occurs in the short-time sector corresponding to relaxation within the bulk 
of domains. In this X(t,s — > t ) — » 7 + n — > 1, the identity 7 = fj, gives 

Xoo = 7 = ji = 1/2, a result that has been found in the one-dimensional Glauber-Ising 
chain [23 El El- Let us mention that FD plots similar to those shown here (Figs. ^ 
EJ) with a limiting value = 1/2 have been reported in one-dimensional models |12j . 
yet it would be interesting to rigorously explain their emergence along the present line of 
reasoning. The conclusion we draw from this paragraph is that a scenario where X^ is 
finite at T = does not appear compatible with an entropically driven relaxation. 

The validity of the expression Wh ~ w exp(— 0hfi) requires that all rates corresponding 
to all elementary activated processes are modified by the same amount (up to linear order 
in h) by the perturbation, a requirement that may not be always fulfilled. For instance, 
let us consider a spin-glass Ising system staying in a locally stable energy minima while 
aging (i.e. the configuration is stable against all one-spin flips). Assume also that P spins 
(pointing in different directions) need to be simultaneously activated for the system to 
escape from that local minima. If there is no long-range order (as is common in spin 
glasses) the local field distribution is symmetric around zero, so the local fields acting on 
these P spins will point in different directions. A uniform perturbing field will therefore 
increase or decrease the local field, and the modified rates Wh ~ wexp(-/?/j/j) will not 
be the same for all P spins. Therefore, the mechanism behind the emergence of effective 
temperatures in frustrated glassy systems (as compared to coarsening systems) seems better 

9 Several scenarios based on this cooperative behavior have been proposed 1181 1191 to explain the origin of 
super-activated effects in structural glasses. 
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described by the BMM rather than by the BM. 

There are two additional aspects we would like to discuss. One aspect concerns the 
predicted waiting-time dependence of Xoo and T e g in IS models. For trap models described 
by a time-independent density of traps QlOf) such dependence is not expected, indeed the 
limiting values (|31I42|) are time independent. For generic glassy systems, though, the 
width of the exponential density of traps (i.e. 7 according to (J2HI) should be time 
dependent as the average energy decreases. Including such dependence would produce FD 
plots identical to those found in other solvable models like the oscillator model The 
second aspect concerns the characteristic form of FD plots in realistic glass models. How 
to reproduce IS FD plots with two linear segments, one for C — > 1, the other for C — ► 0? 
This requires the introduction of some internal dynamics within traps which is absent in 
the original formulation of both the BM and the BMM. In the same vein, generalization of 
this treatment beyond the IS pattern to many timescales could be accomplished through 
more general transition rates depending on trap correlations q a , a ' = M a M a i in some type 
of hierarchical structure @]. 

Finally some words about the violation of the property of detailed balance that we 
have tacitly assumed in the rates I|23I34|> . Detailed balance is of course violated in the 
BMM at T = 0. As we have already insisted on (footnote 0} this assumption is essential, 
as trap models are coarse-grained descriptions in phase space of a microscopic dynamics 
that satisfies that property. The rates W(a'|a) in © are obtained by projection of the 
microscopic probabilities and rates onto the corresponding coarse-grained representation 
of the original phase space. In this projected description detailed balance is not fulfilled 
anymore (for a discussion see Sec. 3. 3 in 6 ). Moreover, detailed balance is known to be 
violated in driven systems, yet they also show IS FD violating patterns. Therefore, an 
explanation on the emergence of effective temperatures cannot rely on the fulfillment of 
the detailed balance property but need to be more general. A coarse-grained description 
of glassy dynamics, where rates are entropically driven, accomplishes this goal. 

6 Concluding remarks 

Trap models represent a useful approach to understand universal aspects of glassy dy- 
namics. They can be seen as a coarse-grained description of structural relaxation where 
nucleation of cooperatively rearranging regions is modeled by jumps among traps. They 
are expected to be relevant to describe relaxation processes in structural glasses. The vast 
majority of theoretical and numerical studies on glassy systems show that one-step (IS) 
FD plots are common in the aging regime. However, recent studies on trap models have 
not found indications of IS behavior in FD plots, casting doubts on the suitability of trap 
models to reproduce this particular feature of glassy systems. 

In this paper we have addressed this issue by showing that proper treatment of the 
influence of the perturbation on the transition rates partially resolves the problem. Indeed, 
IS plots can be obtained if entropic effects, driven by the unbalance of the population 
of traps in the presence of a field, are included in the analysis. We concentrated our 
attention in both the Bouchaud model (BM) and the Barrat-Mezard model (BMM), in 
this last case at zero-temperature to better emphasize the main result. We find that 
Xoo = lim^oo X(t, s) (for the BM) and = lim^^ T cS (t, s) (for the BMM at T = 0) 
are universal quantities that only depend on a parameter (7) describing the dependence 
of the perturbed rates on the observable value at the arrival trap. Yet, the physical origin 
of this dependence appears very different in both models. While perturbed rates are 
activated in the BM, they are not in the BMM but originate from the entropic unbalance 
of traps upon the action of the field. This shows that IS behavior in glassy systems, and 
the emergence of finite effective temperatures in the low T limit, are the fingerprint of 
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a dynamics governed by entropic relaxational processes. Straight IS plots (compatible 
with a unique long-time sector) are obtained in both models only when escape rates are 
unchanged (up to the linear order) by the action of the field. 

Which model (BM or BMM) represents more faithfully the slow dynamics of IS sys- 
tems? The most probable answer is that neither the BM or the BMM provides a complete 
qualitative description of glassy relaxation and that aspects of both models must simul- 
taneously intertwine in the overall relaxation. The dependence of the perturbed rates on 
the observable taken at the departure trap (the term including the parameter /i in (j2HJl) 
can be identified with energy barriers governing the thermally activated escape process 
from a trap. While the dependence on the observable taken at the arrival trap (the term 
including the parameter 7 in H34I) ) is consequence of the unbalance of the number of traps 
induced by the field, leading to temperature independent entropic relaxation. Also, the full 
straight-line FD plots obtained for /1 = represent an extremely simplified behavior of a 
more realistic scenario where a dependence of the trapping time at linear order in the field 
is expected. Yet, this dependence is not in conflict with the emergence of "quasi-straight" 
FD plots for which the existence of a limiting value Too(s) is enough. 

There remain several open question for the future. It could be instructive to extend 
this analysis to the BMM at finite temperature. However, a relevant question is to under- 
stand how all these considerations are to be generalized to include the case of arbitrary 
observables. For the non-neutral case, explicit computations as have been done here are 
not possible. Only for the purely activated case calculations can be done JT], but that 
analysis does not include the entropic effects discussed in this paper. Yet, there may exist 
general proofs from which one can deduce specific predictions to be challenged in several 
microscopic theoretical models and well as in numerical simulations. 

Acknowledgments. I warmly thank P. Sollich for communicating me his results after 
submission of the first draft of the paper and the fruitful exchange of ideas that followed 
after. FR is supported by the Spanish Ministerio de Ciencia y Tecnologia Grant BFM200T 
3525 and Generalitat de Catalunya. FR has benefited from the warming hospitality at the 
Bustamante lab in the University of California in Berkeley where this work has been done. 

Note added: Just after completion of this worked I learned of calculations by Sollich 
showing that, for the BMM, purely activated rates 1)11)) lead to FD plots not displaying 
IS behavior, yet they seem to be described by \34\j with 7 = 0, fi = 2/T g , i.e. without 
dependence of the rates upon the value of the observable taken at the arrival trap. FD plots 
are therefore very similar to those found in the BM with purely activated rates. This result 
is very interesting as it shows that the emergence of IS behavior in trap models cannot be 
described by only shifting the energies according to the value of the observable [cfr. 111)) ]. 
Identification of mechanisms behind the emergence of perturbed rates such as )2'^4\j , 
beyond the purely activated ones, is certainly needed to clarify this issue. 
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